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We have studied the high temperature multiferroic cupric oxide using polarized neutron diffraction 
as a function of temperature and applied electric field. We find that the chiral domain population 
can be varied using an external electric field applied along the b axis. Using representation anal- 
ysis we derive the incommensurate magnetic structure in the multiferroic phase. The origin of 
the magnetoelectric coupling is consistent with models based on the inverse Dzyaloshinskii-Moriya 
interaction, but is different from the simple cycloidal mechanism. 
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I. INTRODUCTION 

Materials in which ferroelectric and magnetic order pa- 
rameters couple are of great interest for application in 
spintronic devices. Such materials are termed multifer- 
roic and have been actively studied in recent years.— ~— Of 
particular interest are the type-II multifcrroics in which 
the magnetic order is directly and strongly coupled to the 
electric polarization. The strong coupling is often associ- 
ated with non-collinear magnetic structures which orig- 
inate from competition and frustration in the magnetic 
interactions. Most of the strongly-coupled multiferroics, 
therefore, exhibit their unusual properties at low tem- 
peratures (< 40 K). An exception is cupric oxide, which 
develops a multiferroic phase below 230 K A 

The magnetic phase diagram of CuO in zero magnetic 
field contains two phases, AF1 and AF2£~— The low tem- 
perature paraelectric phase (AF1) is a simple commensu- 
rate antiferromagnet which is stable up to Tni = 213.7 K. 
At Tni there is a first-order transition to the multiferroic 
phase (AF2), which has a complex, non-collinear incom- 
mensurate magnetic structure and a ferroelectric polar- 
ization along the crystallographic b axis. Both magnetic 
and ferroelectric orders vanish above Tn2 = 230 K. 

The mechanism that causes the simultaneous emer- 
gence of ferroelcctricity and magnetic ordering in CuO is 
still under discussion and a number of theories have re- 
cently been proposed. A phenomenological Landau the- 
ory approach was employed by Toledano et al. (Ref. [l(]) 
to investigate the two ordered phases. Using this theory 
they were able to explain the sequence of phases, with 
the incommensurate multiferroic phase at a higher tem- 
perature than the simple antiferromagnetic phase, and 
predicted that both phase transitions are first order. Gio- 
vannetti et al. (Ref. Illh and Jin et al. (Ref. [l2l ) performed 
density functional calculations which suggest that the fer- 
roelectric polarization is induced via the Dzyaloshinskii- 
Moriya (DM) interactio n 13 ' 14 between neighboring chains 
of non-collinear spins. They further argue that the small 
incommensurability in the AF2 phase does not play an 



important role in the multiferroic mechanism. 

Overall, it is widely accepted that the competing ex- 
change interactions together with geometrical frustration 
lead to a non-collinear magnetic ordering which breaks 
inversion symmetry and allows polar lattice distortions. 
What is different about CuO is that a strong Cu-O-Cu 
superexchange interaction along the [1,0, 1] direction al- 
lows it to remain ordered at high temperatures^^ - — 

If the magnetoelectric coupling is strong enough, it 
should be possible to vary the magnetic domain popu- 
lation with an external electric field. The domain pop- 
ulation can be monitored by polarized neutron or pho- 
ton diffraction. This type of experiment has been per- 
formed previously on several different multiferroics, and 
as well as demonstrating the manipulation of magnetic 
domains by electric fields these experiments have pro- 
vided useful information about the magnetoelectric cou- 
pling mechanisms *£r— 

Here we present the results of a polarized neutron 
diffraction study of the magnetic response of CuO in an 
applied electric field. We demonstrate that the magnetic 
domains in CuO can be switched by an electric field ap- 
plied along the ferroelectric polarization axis. The field 
required to reverse the domain population is found to 
be closely related to the temperature. However, it does 
not appear to be possible to transform the crystal into a 
single magnetic domain even at large applied fields. 



II. EXPERIMENTAL PROCEDURE 

A single crystal sample of CuO was grown from a 
melt of high purity cupric oxide (99.995%) by the opti- 
cal floating-zone method^ The crystal, which was single 
phase and untwinned, was cut with parallel flat faces per- 
pendicular to the b axis giving a plate-like sample of sur- 
face area of 21mm 2 , thickness 0.9mm and mass 0.175 g. 
An 8-10 nm layer of chromium and 40-50 nm of gold was 
evaporated onto the flat faces. Electrical contacts were 
attached using silver paste and the sample was sealed in 
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FIG. 1. (color online). The observed magnetic structures of 
CuO in projection onto the a-c plane. The four Cu sites in 
the conventional structural unit cell (indicated by the paral- 
lelogram) are labelled 1-4. The colored circles represent the 
component of the moment along b: light red is along +b, dark 
blue is along — b, and blue circles with arrows are moments 
which lie in the a-c plane. Panel (a) shows the commensu- 
rate magnetic structure in the AF1 phase whose moments are 
all parallel or antiparallel to b. Panel (b) shows the spin ar- 
rangement in the incommensurate AF2 phase. The moments 
lie in the planes indicated by the dotted lines, which contain 
the 6 axis and the direction v. The AF2 structure has two do- 
mains of opposite handedness. Only the left-handed domain 
is shown here. 

an aluminium vacuum can to avoid electric field break- 
down. 

Polarized neutron diffraction measurements were per- 
formed on the TASP/MuPAD instrument at SINQ, Paul 
Schcrrer Institut^^ The sample was aligned with the 
a and c axes in the horizontal scattering plane such that 
the applied electric field was along the —b direction. The 
scattering vector is defined as Q = k; — kf where k; 
and kf are the incident and scattered neutron wavevec- 
tors, and fcj = fcf = 1.97 A -1 . The flipping ratio, as 
measured on the (002) and (200) structural Bragg peaks, 
was found to be 17. The experimental data presented 
here have not been corrected for the non-ideal neutron 
polarization, but corrections have been included in all 
calculations from magnetic structure models. The back- 
ground signal in the different polarization channels was 
measured at positions away from the Bragg peaks and 
found to be negligible. 

Pyroelectric current measurements to obtain the elec- 
tric polarization of CuO were made with a customized 
PPMS (Quantum Design) insert. The crystal was cut 



from the same rod as the neutron diffraction sample. The 
surface area of this sample was 10 mm 2 and its thick- 
ness was 1.7 mm. The sample was cooled from 250 K in 
the paraelectric phase to 220 K in the multiferroic phase 
with a potential difference of 180 V applied parallel to 
b. Upon reaching 220 K, the electric field was turned off 
and the surface charge on the Cr/Au coating was allowed 
to discharge. The pyroelectric current was then recorded 
either upon warming to above Tn2 or cooling down be- 
low Tni , using the same field-cool procedure before each 
measurement. The heating/cooling rate was lKmin -1 . 
To eliminate systematic errors, an applied electric field 
parallel and antiparallel to b was employed. 

Zero-ficld-coolcd magnetization measurements were 
performed with a SQUID-based magnetometer (Quan- 
tum Design). The measuring field was 1000 Oe. A slow 
heating rate of 0.01 KmmT 1 was used close to the sharp 
Tni transition, and a rate of 0.1 K min -1 was used above 
and below it. 



III. MAGNETIC REPRESENTATION 
ANALYSIS 

The antiferromagnetic structure of the AF1 phase 
which exists below Tni ~ 213 K is shown in Fig. [Ha). 
The magnetic propagation wavevector is q cm = 
(0.5, 0, —0.5), and the spins lie parallel to the b axis. The 
multiferroic phase (AF2) exists within the narrow tem- 
perature range between Tni and Tn2 ■ Neutron diffraction 
experiments have shown that the magnetic structure in 
this phase has an incommensurate propagation wavevec- 
tor q icm = (0.506, 0, -0.483). The Cu 2+ moments rotate 
with an approximately circular envelope in the plane con- 
taining the vectors b and v = 0.506a* + 1.517c*, where 
a* and c* are reciprocal lattice basis vectors*^ Spon- 
taneous ferroelectric polarization is found to develop in 
the AF2 phase along the b axis^ No ferroelectricity is 
observed for T < T N1 or T > T N2 . 

The crystal structure of CuO can be described by the 
monoclinic C2/c space group, with /3 = 99.55° and lat- 
tice parameters of a = 4.69 A, b = 3.42 A and c = 5.13 A 
at 300 Kj2! For both magnetic phases, the magnetic 
propagation vector q is left unchanged by the identity 
,9i = {E\0, 0,0} and glide plane <?2 = { cr y 1 , 0, 0.5} sym- 
metry operations of the paramagnetic space group. The 
little group G q is therefore {(71,32}- The irreducible rep- 
resentations of G q are Ti and T2, whose characters are 
given in Table HI The primitive cell of CuO contains two 
Cu atoms, but the G-centering {T|0.5, 0.5, 0} then gener- 
ates two more Cu sites such that there are four Cu atoms 
in the conventional unit cell. 

The components of the two Cu spins in the prim- 
itive cell form a 6-dimensional representation r mag of 
G q which reduces to r mag = 3Fi + 3r2. Symme- 
try analysis^ fixes the relative phases between Cu sites 
and the experimentally-determined spin directions fur- 
ther constrain the magnetic structures. The magnetic 
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TABLE I. Character table of the irreducible representations 
of G q with the phase factor given as tp = — exp(inq z ). 
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M b 



Mr 



1 cos(q • t + 9b) cos(q-t + #!,) 

2 cos(q ■ t + q c TY + 6 b ) - cos(q • t + q c 7r + 9 V ) 

1 cos(q • t + 9b) cos(q ■ t + #„) 

2 - cos(q ■ t + q c TV + 6 b ) cos(q ■ t + q c n + 9 V ) 



TABLE II. The symmetry-allowed magnetic structures 
of CuO. The AF2 structure has ordering vector q = 
(0.506, 0, —0.483) and can be resolved into components par- 
allel to b and v = 0.506a* + 1.517c*. The AF1 structure 
has q = (0.5, 0, —0.5) and only has components along the 6 
axis. The Cu site positions j are: (1) |,j,0; (2) |,|,|; (3) 
§,§,0; (4) \,\,\, as indicated in Fig. [T] The lattice trans- 
lation vectors are denoted by t, and the moments of atoms 
on sites 3 and 4 are found using the C-centering translation 
tc = (1/2, 1/2, 0), resulting in an additional phase of q • tc- 
For the AF1 structure, 9b = tt/4, and for the AF2 structure, 
9 b = and 9 V = tt/2. 



different symmetry for the b and v magnetic components 
(A/ b (1) Mi 2) or M^MP) do have a net electric polariza- 
tion in the spin current model. 

Of the four assumed AF2 magnetic structures for CuO, 
only MfrMv is consistent with neutron diffraction ex- 
periments. The structure can be ruled out be- 
cause it predicts a very small structure factor (identically 
zero for an isotropic Cu form factor) for the Q = qi cm 
fundamental diffraction peak, while the M^p My 1 ^ and 

(2) (2) 

M b My structures give very poor agreement with neu- 
tron polarimetry measurements^ The M b My struc- 
ture is shown in Fig. [T^b), arL d has been used in all sub- 
sequent analysis presented here. 

By enumerating the possible symmetry-allowed mag- 
netic structures for the AF1 and AF2 phases of CuO and 
comparing these with experiment, we conclude that the 
structures shown in Fig. [T] are the correct ones. We can 
therefore use these to quantify changes in domain pop- 
ulation induced by an applied electric field. In the next 
section we describe briefly the methodology of our exper- 
iment. 



IV. SPHERICAL NEUTRON POLARIMETRY 



components of the structures consistent with this analy- 
sis are given in Table [TTJ 

The commensurate AF1 phase can be described us- 
ing just M&, the basis vector along b. The spin struc- 
tures corresponding to the pure Ti and T2 irreducible 

representations, Mi and M b , are given in Table ILTl 
The choice of phase factor Oh = 7r/4 fixes the amplitude 
of the magnetic moment on each Cu site to be equal. 7 
Single-crystal neutron diffraction measurements of mag- 
netic Bragg peak intensities below Tni have established 
that the structure for the AF1 phase has the T\ symme- 
try. 

The incommensurate AF2 phase has magnetic compo- 
nents along both the b and v directions, M b and M v , re- 
spectively. To correctly describe the AF2 spin structure 
we set 9h = and 9 V = tt/2 so that the spins rotate with 
a circular envelope, as observed experimentally. Assum- 
ing the Mb and M v components can each be described 
by a single irreducible representation, Ti or T2, there 
are four possible magnetic structures which we denote 
by M b {n) M v m) . These structures are chiral, so each has 
two domains of opposite handedness related by spatial 
inversion. When the b and v magnetic components be- 
long to the same irreducible representation, i.e. M^ M v ^ 
(2) (2) 

or Ml M v , this corresponds to a spin arrangement in 
which the spins on sites 1 and 3 rotate in the opposite 
sense along the a axis to those in sites 2 and 4. This 
results in zero net electric polarization in the unit cell 
assuming the magnetoelectric coupling depends on the 
spin current Si x S2. The two magnetic structures with 



In general, magnetic scattering of neutrons involves a 
change in the polarization state of the neutrons. Spher- 
ical neutron polarimetry is a technique which measures 
the polarization of the scattered neutron beam for an 
arbitrary incident beam polarization. The results are ex- 
pressed in terms of a polarization matrix, 
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where a and /3 indicate the polarization of the incident 
and final beam, respectively. I a p is the magnetic Bragg 
peak intensity, given by 

I a p = |(a|Mj_(Q) • <t|/3)| 2 , (2) 

where er is the Pauli operator for the neutron spin and 
Mj_(Q) is the component of M(Q) perpendicular to Q, 
where M(Q) is the Fourier transform of the magnetiza- 
tion. 

For CuO, M(Q) may be written, 



M(Q) = EEww( t ) eiQ ' 



(t+r 



(3) 



where rj is the displacement of the jth Cu site (j = 1-4) 
from the origin of the conventional unit cell, /j(Q) is the 
Cu magnetic form factor, and t is a lattice translation 
vector. The components of Mj-(t) are given in Table HT1 
In an ideal measurement with perfect neutron polar- 
ization the polarization matrix for CuO is given by 
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where, 
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Here we have used the Blume coordinate system , 30 ' 31 
with the x axis parallel to Q, the z axis perpendicu- 
lar to the scattering plane, and the y axis chosen to 
make up a right-handed set of axes. In this axis sys- 
tem = (0, M y ,M z ). The existence of a non-zero 
'chiral' term C in the polarization matrix implies a non- 
collinear magnetic structure that is, by definition, con- 
structed from orthogonal, out-of-phase components, such 
as a helix or cycloid. Therefore, C = in the collinear 
AFl structure of CuO, and in the AF2 phase, left- and 
right-handed chiral domains will give C < and C > 0, 
respectively. Equivalently, if we measure the P yx and P zx 
elements on a magnetic reflection we can determine the 
relative proportions of the two chiral domains. 

The parameters C and D can be calculated from Eqs.[5] 
for a general helicoidal magnetic structure. For the par- 
ticular case when the spins rotate with a circular envelope 
in a plane perpendicular to q, the values for the reflec- 
tion Q = q arc C = ±1 and D = 0. In the AF2 phase of 
CuO, the propagation vector q; cm is not perpendicular 
to the plane of rotation of the spins, but rather lies at an 
angle of 107° to it — Fig. [T^b). For Q = q this results 
in a small but non-zero value of D, and a magnitude of 
C which is slightly less than 1 for a single chiral domain. 

The polarization matrix measured on a nuclear reflec- 
tion contains only the diagonal elements P Q( g = 6 a p. 
Measurements on nuclear (or equivalently magnetic) re- 
flections are useful in estimating the amount of depolar- 
ization of the incident neutron beam due to stray fields. 
This is quantified by a flipping ratio, defined for a nu- 
clear Bragg peak as the ratio of the scattering inten- 
sity measured in the non-spin-flip to spin-flip channel, 

R I eta. / Icta • 



V. EXPERIMENTAL RESULTS 

Temperature profiles of magnetic and dielectric prop- 
erties of our CuO crystal are presented in Fig. [2] Fig- 
ure [5Ja) shows temperature sweeps of the magnetic sus- 
ceptibility recorded along the reciprocal lattice axes a*, 
b* and c* . A sharp discontinuous transition is observed 
along all directions at Tni = 213.7 K. Above Tni the 
susceptibility increases linearly with temperature up to 
Tn2 ~ 230 K at which point a change in slope is observed. 
Our data are in very good agreement with previous mea- 
surements on CuO A Susceptibility measurements were 
made on cooling and warming. We find no hysteresis 
at Tni to within 0.2 K, despite the first-order nature of 
this transition. We observed a similarly abrupt and non- 
hysteretic transition at Tni in magnetic neutron diffrac- 
tion measurements also made on cooling and warming. 
This is at odds with the large hysteresis of approximately 
5 K reported in diffraction measurements by Yang et al. 
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FIG. 2. (color online). Magnetic and dielectric properties 
of CuO. The upper panel (a) The zero-field-cooled magnetic 
susceptibility measured with a field of H = 1000 Oe applied 
parallel to the {a*,b*,c*} reciprocal lattice vectors, (b) 
The electric polarization with the electric field applied along 
b* direction. Measurements were made with a poling field 
of 106 V mm -1 applied in opposite directions (red and blue 
curves), (c) shows the temperature dependence of the coer- 
cive field measured by neutron polarimetry. 



(Ref.i). The different behaviors may relate to slight dif- 
ferences in the chemical composition of the CuO crystals. 

The electric polarization P e along ±b* obtained by in- 
tegrating the pyroelectric current with a poling field of 
106 V mm -1 is shown in Fig. [UJb). A sharp depolariza- 
tion is observed close to 213 K coincident with the AF1- 
AF2 magnetic transition. A maximum polarization of 
approximately 100 /zC m -2 is recorded. 

As discussed above, neutron polarimetry can be used 
to determine the relative proportions of the two chiral 
domains in CuO. Any changes in population induced by 
an applied electric field will affect the P yx and P zx ele- 
ments in the polarization matrix while leaving the other 
elements unchanged. 

Complete polarization matrices P a p measured with an 
applied external electric field of ±670 V mm -1 , recorded 
at 220 K, are shown in Figs. and (b). The data 
clearly shows a reversal in the sign of the terms P yx 
and P zx . A decrease in the xx component in Fig. [3jc) 
is anomalous. Its origin is unclear but may be due to 
a small misalignment of the crystal with respect to the 
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to be non-zero. 

Measurements of the P yx component of the polariza- 
tion matrix at a constant temperature of 220 K after ini- 
tial field cooling are shown in Fig. Etc). We find that 
by sweeping the electric field, which is applied in the - 
b*dircction, from670Vmm _1 to — 670Vmm _1 and back 
to 670 Vmm -1 we observe a change in sign and magni- 
tude of P yx . It follows that the electric field must be cou- 
pled to the chiral magnetic domains. The coercive field, 
defined as the electric field at which the chiral domains 
are equally populated, is approximately 90Vmm _1 at 
this temperature. 

We measured P yx loops as a function of electric field at 
a series of temperatures within the AF2 phase. The co- 
ercive field obtained from these loops, shown in Fig[5Jc), 
are found to increase with decreasing temperature. The 
system becomes softer at higher temperatures and there- 
fore the electric field required to balance the domain pop- 
ulation is reduced up to the point of the phase transi- 
tion at 230 K. No hysteresis and hence no coercive field 
is found in the AF1 or paramagnetic (PM) phases. To 
verify that the electric field switching of domains exists 
only in the AF2 phase, hysteresis loops were also made 
at 240 K (centered on qi cm ) and 210 K (centered on qj cm 
and q cm ). These measurements showed field-independent 
behavior of P yx , as expected. 

We note that the magnitude of P yx never reaches the 
maximal value of close to ±1, as expected for the case 
of a single chiral domain in Eq. [5] Even at large electric 
fields a significant depolarization is observed such that 
a saturation value of \P yx \ ~ 0.5, is reached at around 
±500 V mm -1 when measured at 220 K. Hysteresis loops 
measured at both lower and higher temperatures within 
the AF2 phase show similar behavior with the saturation 
value of \P yx \ ~ 0.5. 

To examine this effect further we performed rocking 
scans through the qi cm Bragg peak. Figures E^d) and 
(e) show the intensities in the spin-flip and non-spin-flip 
yx polarization channels at the extrema of the hysteresis 
loop. We would expect the left-handed helical domain 
(i.e. that shown in Fig. [TJb)) to scatter mainly into the 
yx polarization channel, and the right-handed helical do- 
main to scatter mainly into the yx channel. Taking into 
account the non-ideal neutron beam polarization in the 
experiment we calculate the ratio I yx /I yx to be 15 for 
a single magnetic domain. However, from Fig. |3]Jd) this 
ratio is closer to 3.4. 

This discrepancy suggests either that the assumed 
magnetic structure for the AF2 phase of CuO is incorrect, 
or that complete switching between magnetic domains is 
not achieved even though P yx apparently saturates at the 
highest measured electric fields. The magnetic structure 
in the AF2 phase of CuO has been tightly constrained by 
previous neutron diffraction and polarimetry studies^^ 
as well as analysis of our own results,— and so it seems 
reasonable to assume that the magnetic structure that 
we consider is indeed the correct one. Moreover, it can 
be shown that P yx + P yy + P yz = 1 for a single domain of 
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FIG. 3. (color online). Spherical neutron polarimetry data 
for CuO with an applied electric field, (a) and (b) represent 
the complete polarization matrix at qi cm = (0.506, 0, —0.483) 
with applied field E of — 670 Vmm -1 and +670 Vmm -1 , re- 
spectively. The color bars are the measurements, and the 
symbols (o) are calculated assuming theM 6 (1) Mi 2) AF2 mag- 
netic structure and taking into account the non-ideal neutron 
beam polarization. The fitted populations of the two chiral 
magnetic domains are indicated. Panel (c) shows the P yx 
component of the polarization matrix during an electric field 
sweep at a fixed temperature of 220 K. Panels (d) and (e) show 
rocking scans through the qi cm magnetic Bragg peak as mea- 
sured in the yx and yx polarization channels. Measurements 
were made at constant applied electric field of —670 V mm -1 
(panel d) and +670 Vmm -1 (panel e). 



any magnetic structure. This sum rule is not satisfied by 
our polarimetry data — see Figs. EJa) and (b). We con- 
clude, therefore, that electric field conversion to a single 
magnetic domain is never achieved in CuO. 

To quantify the domain populations, we have used the 
model for the AF2 magnetic structure where basis vectors 
along b correspond to the ri representation and the basis 
vectors in a-c plane to T2. A correction was made to the 
simulations to account for non-ideal beam polarization 
assuming a beam polarization efficiency of 94%. We fit- 
ted the complete measured polarization matrices allow- 
ing just the domain fraction to vary. The domains Dl 
and D2 correspond to left- and right-handed spin struc- 
tures, respectively. The electric field switches between 
the domains, with a large positive electric field yielding 
a predominantly left-handed domain and, conversely, a 
large negative field promoting the right-handed domain. 

We find that typically, in the AF2 phase, at electric 
fields approaching saturation in P yx , the magnetic do- 
mains arc populated in approximately 80:20 proportion. 
We have observed, therefore, that the electric field does 



6 



not induce a single magnetic domain. Possible reasons 
for this are, (i) crystal defects acting to pin the magnetic 
structure and prevent full domain reversal, (ii) the high- 
temperature of the multiferroic phase causing thermal 
relaxation of the domain population in the time frame of 
the experiment, or (iii) a reduced electric field strength 
near the edges of the crystal caused by the incomplete 
coating of the surfaces with the gold electrodes. 

VI. DISCUSSION 

Density functional theory calculations of the electric 
polarization in the multiferroic phase of CuO were re- 
ported in Refs.[rHandll2l. The magnitude of P e predicted 
by both groups is consistent with the experimental bulk 
value of approximately 100 /^C m~ 2 . However, we observe 
that the direction of P is along +b for the left-handed 
domain, in contradiction to the direction predicted by 
one of the models*^ 

The magnetically-induced ferroelectric polarization in 
CuO has been explained through a number of models, 
all of which are based upon the inverse Dzyaloshinkii- 
Moriya interaction between nearest-neighbor, non- 
collinear spins. For such spin structures, the DM inter- 
action may reduce the system's energy, as described by 
the following term in the Hamiltonian, 

7/,,m ^D, ;S • S, :. (6) 

The DM vector, Dy, can be resolved into components 
perpendicular and parallel to the vector r,j connecting 
the spins as, 

T)ij = Pjj x Vij + (JijYij, (7) 

where Py and aij arc a polar vector and pseudoscalar, re- 
spectively. In CuO, Pij is proportional to a local electric 
polarization, and is related to the magnetic chirality. 
Usually only the first term is considered when explaining 
the evolution of ferroelectricity, as it explicitly includes 
a polar vector. This is indeed the case in all previous 
studies of CuO. Kimura et al. (Ref. consider a spin 
cycloid whose propagation vector is along the direction 
of the incommensurate modulation qi cm — Qc.™, and Gio- 
vanetti et al. (Ref. ITT1 ) and Jin et al. (Ref. fl3 ) find that 
the important magnetoelectric interaction is between ap- 
proximately perpendicular spins on neighboring chains of 
Cu atoms running along the [1,0,1] direction and does 
not depend on the incommensurate modulation. 

The second term of Eq. [7] does not inherently give rise 
to a polarization. However, it has recently been shown 
that the chirality of a magnetic structure, a^, can in- 
duce an electric polarization through coupling to a unique 
structural rotation^ Furthermore, the electric polariza- 
tion is then strictly constrained to lie parallel to the rota- 



tion axis. This coupling is limited to a small class of crys- 
tal structures which are termed ferroaxial.— CuO adopts 
the space group C2 /c, a member of the ferroaxial crystal 
class, with a two-fold rotation axis parallel to the crys- 
tallographic 6-axis. Coupling between magnetic chiral- 
ity and the crystal structure could therefore result in an 
electric polarization parallel to the b observed. It 

would be of interest to include this chiral term in future 
DFT calculations. 

A recent study of the bulk magnetoelectric properties 
of polycrystalline CuO showed that the application of 
a magnetic field has little effect on the ferroelectricity 
and concluded that the magnetoelectric coupling is very 
weaki^ Our method of applying an electric field and 
measuring the effect on the magnetic structure shows 
that the magnetoelectric coupling is in fact strong. The 
difference between these experiments is that in our work 
we are probing the coupling between the macroscopic po- 
larization and a spatially- varying magnetization, whereas 
in Ref. [H the coupling is between the macroscopic polar- 
ization and a uniform magnetization. These studies are 
therefore complementary, not contradictory. 



VII. SUMMARY 

We have used polarized neutron diffraction to show 
that an electric field applied along the b axis is able 
to switch between magnetic domains in CuO. The re- 
sults demonstrate that the ferroelectric polarization is 
directly coupled to the chiral magnetic order. We have 
solved the magnetic structure in the multiferroic phase 
by representation analysis and confirmed that it is con- 
sistent with previous data and that it supports a ferro- 
electric polarization. The microscopic origin of the ferro- 
electricity is consistent with models based on the inverse 
Dzyaloshinkii-Moriya mechanism^ - — 
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